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STATEMENT  OF  PROBLEM 


Printed-circuit  or  microstrip  antennas  have  proven  to  be  extremely 
useful  in  a  wide  variety  of  applications  requiring  radiators  that  are 
both  low-profile  and  conformal.  Previous  theoretical  Investigations 
have  relied  on  simple  models,  which  provide  very  useful  information, 
but  are  limited  in  their  applicability  and  accuracy.  During  the 
present  project  a  more  rigorous  theoretical  approach  has  been  developed. 
In  addition  a  radiator  with  a  new  geometrical  shape,  that  of  an  ellipse, 
has  been  investigated  both  theoretically  and  experimentally.  The 
resulting  antenna  is  able  to  provide  circular  polarization  using  only 
a  single  simple  feed  network.  _ _ 


SUMMARY  OF  RESULTS 


Several  problems  associated  with  low-profile  antennas  were 
addressed  in  this  project.  A  large  portion  of  the  effort  was  directed 
toward  printed-circuit  or  microstrip,  type  antennas. 

From  the  time  of  their  inception  printed-circuit  antennas  have 
been  plagued  by  a  seemingly  insurmountable  problem  of  excessively  small 
bandwidth.  From  all  indications  there  seems  to  be  little  one  can  do 
to  improve  the  bandwidth  beyond  simply  increasing  the  thickness  of 
the  substrate.  This  method  cannot,  however,  be  extended  too  far 
without  the  loss  of  the  highly  desirable  low-profile  characteristic 
of  the  antenna.  In  some  applications  increased  bandwidth  is  only 
needed  so  the  antenna  can  operate  at  two  distinct  frequencies  which 
may  be  too  far  apart  for  one  antenna  to  operate  efficiently  at  both 
frequencies.  For  this  purpose  a  dual  frequency  antenna  was  designed 
by  stacking  two  slightly  different  sized  circular  discs.  Through  a 
proper  choice  of  the  two  disc  diameters  and  their  spacings,  the  two 
resonant  frequencies  can  be  adjusted  to  the  desired  values.  The 
technical  details  have  been  previously  published  and  are  included  in 
Appendix  A. 

Early  analysis  techniques  for  printed-circuit  type  antennas 
treated  the  antenna  as  a  resonant  cavity  bounded  above  and  below  by 
conducting  plates  and  on  its  sides  by  perfectly  conducting  magnetic 
walls.  This  assumption  limits  the  model  to  very  thin  dielectric  sub¬ 
strates  and  results  in  an  infinite  input  impedance  at  resonance. 

During  this  project  the  antenna  was  modeled  as  a  cavity  with  a  finite 
admittance  wall  on  its  side  surfaces.  This  admittance  surface  was 
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determined  by  the  radiated  power  and  fringing  field  at  the  edge  of  the 
radiator.  Using  this  improved  model  the  input  impedance  can  be  cal¬ 
culated  directly  and  results  in  better  agreement  with  previous  experi¬ 
mental  data.  The  technical  details  and  specific  comparisons  with 
experimental  results  have  been  published  and  are  shown  in  Appendices  B 
and  C. 

In  its  usual  simple  configuration  whether  rectangular  or  circular, 
the  printed- circuit  antenna  produces  a  linearly  polarized  radiation 
pattern.  In  many  applications,  however,  circular  polarization  is 
required.  Through  the  use  of  multiple  feeds  and  phase  shifters  cir¬ 
cular  polarization  can  be  provided,  but  with  the  obvious  increase  in 
the  complexity  of  the  total  radiating  structure.  In  an  attempt  to 
solve  this  problem,  the  elliptical  printed-circuit  antenna  was  inves¬ 
tigated  theoretically.  The  results  show  that  circular  polarization 
can  be  provided  with  a  single,  simple  feed  by  simply  making  the  shape 
of  the  radiator  slightly  elliptical. 

The  boundary  value  problem  to  be  solved  consisted  of  an  elliptical 
cavity  surrounded  by  a  surface  admittance  wall  as  previously  derived 
for  the  circular  disc  case.  Standard  elliptical  coordinates  can  be 
used  and  the  resulting  fields  expressed  in  terms  of  Mathieu  functions. 
To  provide  circular  polarization  the  ellipse  must  have  a  very  small 
eccentricity  and  is,  therefore,  almost  circular.  For  these  small 
eccentricities  approximate  formulas  are  available  for  Mathieu  functions 
in  terms  of  more  frequently  tabulated  Bessel  functions.  Using  these 
calculated  interior  fields  the  total  radiated  fields  can  then  be 


determined.  For  a  proper  choice  of  eccentricity  good  circular  polar- 


lzation  is  predicted  for  a  feed  located  along  a  radius  at  a  45°  angle 
from  the  major  axis  of  the  ellipse.  The  opposite  handedness  polariza¬ 
tion  can  be  generated  by  placing  the  feed  at  -45°  from  the  major  axis. 
The  detailed  theoretical  derivation  has  been  accepted  for  publication 
and  is  Included  in  Appendix  D. 

A  systematic  experimental  investigation  was  also  undertaken  to 
study  the  elliptical  printed-circuit  antenna.  Several  sets  of  ellipses 
were  etched  on  printed  circuit  boards  of  varying  thicknesses.  The 
slightly  elliptical  discs  varied  from  a  minor  to  major  axis  ratio  of 
1.0  (that  of  a  circle)  to  0.96.  The  field  patterns  were  measured  for 
each  disc  at  several  frequencies  near  the  frequency  where  best  circular 
polarization  was  obtained.  In  addition  the  axial  ratio  was  measured 
for  each  antenna  as  a  function  of  frequency.  Finally,  the  input 
impedance  of  each  radiator  was  measured.  The  data  was  compared  to 
available  theoretical  predictions  and  overall  agreement  was  found. 
Detailed  measurements  are  illustrated  in  a  paper  accepted  for  publica¬ 
tion  that  is  included  in  Appendix  £. 
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A  Dual-Frequency  Stacked  Circular-Disc  Antenna 

STUART  A  LONC.  Ml  Mill  It.  II  t  K.  AND 
MARK  I).  WALTON 


Abslidi  l  The  dual-frcqiieney  behavior  of  a  pair  of  stacked  circular- 
disc  printed -circuit  antennas  is  investigated  experimentally.  The  inpul 
impedance  is  measured  as  a  function  of  the  si/es  of  the  discs  with  em¬ 
phasis  placed  on  the  values  of  the  two  resultant  resonant  frequencies. 


I.  INTRODUCTION 

In  recent  years  the  use  and  interest  in  printed-circuit  (or 
micros! rip)  antennas  has  become  widespread.  A  wide  variety 
of  shapes  and  configurations  have  been  developed  into  useful 
radiating  systems.  In  addition  to  being  low  profile,  their  posi¬ 
tive  attributes  include  ruggedness,  ease  of  construction,  and 
low  cost.  Perhaps  the  sole  major  problem  of  the  entire  class  of 
antennas  is  their  seemingly  inherent  narrow  bandwidth.  Little 
progress  has  been  reported  in  overcoming  this  major  obstacle. 
The  most  usual  method  of  increasing  the  bandwidth  is  simply 
to  increase  the  thickness  of  the  dielectric  between  the  radiator 
and  the  ground  plane.  This  method  cannot,  however,  be  ex¬ 
tended  too  far  without  the  loss  of  the  highly  desirable  low- 
profile  characteristics  of  the  antenna.  Such  a  technique  results 
in  only  minimal  increases  in  the  bandwidth  with  the  resulting 
thicker  antenna  still  being  quite  narrow-band  in  the  usual  sense 
(less  than  2  percent).  Some  significant  increases  in  bandwidth 
have  been  reported  for  a  more  complicated  feed  arrangement 
using  two  probes  and  a  balun  1 1 ) . 

For  some  uses,  however,  the  increased  bandwidth  is  actu¬ 
ally  needed  for  only  two  distinct  frequencies  which  may  be 
too  far  apart  for  a  single  antenna  to  operate  efficiently  at  both 
frequencies.  The  behavior  of  the  antenna  characteristics  for 
the  range  of  intermediate  frequencies  may  be  of  little  or  no 
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l  ip  I  Suckril  circular-disc  antenna  ctiuctuiv 


l  ip  2  Smith  that!  impedance  nl  stacked  structure. 


concern,  For  this  purpose  a  dual-frequency  antenna  can  be  de¬ 
signed  by  stacking  two  slightly  different  sized  circular  discs. 
Through  a  proper  choice  of  the  two  disc  diameters  and  their 
spacings,  two  separate  resonances  can  be  found  and  adjusted 
to  the  desired  values  (21 . 

II.  EXPERIMENTAL  MEASUREMENTS 

The  proper  feeding  method  is  not  necessarily  apparent  for 
the  stacked  circular-disc  structure.  To  provide  the  desired  dual- 
frequency  behavior  the  inner  conductor  of  the  coaxial  feed  is 
allowed  to  pass  through  a  clearance  hole  in  the  lower  disc  and 
is  then  electrically  connected  to  the  upper  disc,  as  shown  in 
Eig.  1.  The  actual  antenna  measured  consisted  of  two  circular 
discs  which  were  photoetched  on  separate  microwave  printed- 
circuit  boards  and  then  carefully  aligned  so  that  their  centers 
were  along  the  same  line  perpendicular  to  the  ground  plane. 
The  driving-point  impedance  was  then  measured  as  a  function 
of  frequency  using  a  network  analyz.cr  for  the  slacked  struc¬ 
ture.  The  sizes  of  the  two  discs  and  their  respective  spacings 
-ould  be  varied  and  the  resulting  behavior  of  the  antenna  char¬ 
acteristics  then  noted. 

The  majority  of  the  measurements  were  carried  out  for  a 
.tructurc  with  the  lower-disc  diameter  2a 2  3.78  cm  and  with 

he  two  spacings  7,  =  72  1  0.075  cm  (approximately  1/32  in 
each)  The  dielectric  material  was  teflon-fiherglass  with  a  rela¬ 
tive  permittivity  of  er  --  2.47.  The  upper-disc  diameter  2 a, 
vas  then  varied  in  several  increments  from  3. 5-4.0  cm  with 
he  feed  position  located  near  a  point  one-half  the  radius  from 
the  center  of  the  disc. 

A  typical  impedance,  that  for  2 a,  3.78  cm  (the  same  size 

s  the  lower  disci,  is  shown  in  the  Smith  chart  of  Fig.  2.  In¬ 


tel 

Tig.  3.  Real  and  imaginary  parts  of  impedance  of  slacked  structure 
(a)  2aj  -  3.70  cm.  (b)  7a\  -  3.78  cm.  (c)  2 a,  *  3.85  cm. 

stead  of  the  usual  circle  that  would  be  measured  for  a  single 
circular-disc  printed-circuit  antenna,  two  very  distinct  circles 
are  seen.  Each  of  these  circles  crossing  the  right-hand  real  axis 
corresponds  to  one  of  the  two  desired  resonances.  Riots  of  the 
real  and  imaginary  parts  of  the  impedance  versus  frequency  il¬ 
lustrate  this  dual-frequency  behavior  more  clearly.  Three  typi¬ 
cal  cases  corresponding  to  values  of  2a,  =  3.70,  3.78,  and  3.S5 
cm  are  shown  in  Fig.  3.  Separations  of  10  percent  between  the 
two  resonances  are  clearly  evident.  This  is  far  too  wide  a  separ¬ 
ation  to  be  covered  by  a  single  disc  with  a  simple  coaxial  feed. 
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I  ip  5.  Radiation  patterns  ol  stacked  circular-disc  antenna,  (a)  /  - 
2  83  (ill/.,  (b)/  =  3.1  (.11/  (c)/=  2.9  GHz. 


The  Innclional  behavior  of  these  resonant  frequencies  can 
he  investigated  lor  varying  sizes  of  the  upper  disc.  The  two 
dashed  curves  in  Fig  4  show  t he  upper  (/,,)and  lower  (Jt  ) 
resonant  frequencies  as  the  diameter  of  the  upper  disc  is 
varied.  For  comparison,  the  theoretical  zeroth-order  resonant 
frequency  (Jt))  for  a  single  disc  of  diameter  2a1  is  also  shown 
in  f  ig  4.  In  addition,  the  theoretical  first-order  resonant  fre¬ 
quency  l/j )  which  takes  into  account  the  effect  of  the  finite- 
size  thickness  of  the  dielectric  |3|  is  also  shown  as  a  function 
of  ns  diameter  for  two  thicknesses,  r/j  and  d j  +  d2.  The  lower 
resonance  I//,)  is  relatively  constant,  remaining  near  the  value 
of  a  single  disc  with  la  -  3.78  cm  and  d  -  .075  cm.  The  upper 
resonance  {),,)  is  highly  dependent  on  the  size  of  the  upper 
disc  with  a  functional  behavior  similar  to  a  single,  but  slightly 
larger,  disc  A  summary  of  the  data  taken,  including  those 
cases  shown  in  Fig.  .3,  is  presented  in  tabular  form  in  lable  I 
It  is  thus  possible  to  efiectively  set  the  approximate  value 
of  /,  by  the  choice  of  the  si/e  of  the  lower  disc,  and  then, 
relatively  independently  set  fv  by  a  proper  choice  of  the 
upper-disc  diameter.  Variations  in  the  positions  of  the  reso¬ 


nances  also  result  when  </]  and  d2  are  changed.  A  limited  set 
of  data  for  this  case  is  included  in  Table  1  for  reference 

The  far-field  radiation  patterns  were  measured  in  an  ane- 
choic  chamber  near  the  measured  resonant  frequencies  for  the 
particular  structure  with  2at  =  2 a2  =  3.78  cm  and  d \  =  d2  = 
0.075  cm.  The  results  are  shown  in  Fig.  5  for  each  of  the 
resonant  frequencies  and  for  one  intermediate  frequency.  In  each 
case  the  fields  in  the  two  principal  planes  are  shown,  first  Eg 
for  0  =  0°  (/  -plane)  and  then  Eg,  for  =  00°  (//-plane).  (The 
feed  is  located  at  tp  =  0.)  At  /t  =  2.83  C.Hz  the  fields  are  seen 
to  he  almost  exactly  the  same  as  those  of  a  single  disc  [4| .  In 
each  plane  a  rather  broad  pattern  is  found  with  the  value  of  Eg 
still  finite  in  the  0  -  ')0°  plane  but  with  a  natural  null  for  Eg, 
there.  At  lv  -  3.1  till/  the  haste  pattern  is  still  the  same,  hut 
some  pattern  deterioration  has  begun  to  show  in  the  form  of  a 
small  dtp  at  broadside.  An  additional  measurement  at  /  =  2.9 
(3)1/  shows  that  Ihc  pattern  is  still  very  well  lormed.but  there 
would  be  a  very  sizable  impedance  mismatch  at  this  frequency. 
Il  is  Iheiefore  seen  that  the  narrow  bandwidth  is  primarily  due  to 
the  impedance  variation  and  not  due  to  pattern  deterioration. 
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Pig.  6.  Real  pari  of  impedance  witn  lower  disc  driven. 


III.  STACKED  STRUCTURE  WITH  LOWER  DISC  DRIVEN 

Connecting  the  feed  directly  to  the  lower  disc  and  allowing 
the  upper  one  to  he  purely  parasitic  produces  an  entirely  dif¬ 
ferent  impedance  behavior.  When  plotted  on  a  Smith  chait  the 
impedance  forms  only  one  circle  in  the  same  fashion  as  a 
single  disc.  For  this  radiator  the  inclusion  of  Ihe  upper  disc 
affects  hoth  Ihe  resonant  frequency  and  Ihe  maximum  value 
of  the  resistance.  This  behavior  is  seen  in  Fig.  6  where  Ihe  real 
pail  of  the  impedance  is  plotted  versus  frequency  for  several 
different  sizes  of  upper  discs.  Each  are  edge  driven  with  2a2 
3.7K  cm,  d|  ifj"-0.07S  cm,  and  the  particular  upper-disc 
diameter  as  shown.  I  he  inclusion  of  larger  discs  is  seen  to  shift 
the  resonant  frequency  down  in  value  and  also  to  decrease  the 
value  of  the  maximum  resistance.  Some  slight  improvement  in 
bandwidth  is  also  seen  from  the  relatively  broader  resonance 
curves. 


IV.  CONCLUSION 

When  additional  bandwidth  is  only  needed  at  two  discrete 
frequencies  the  stacked  circular-disc  antenna  can  be  useful. 
Even  though  the  bandwidth  around  each  resonant  frequency 
is  only  about  1  percent,  two  frequencies  separated  by  10  per¬ 
cent  or  more  can  be  accommodated.  A  proper  feeding  struc¬ 
ture  is,  however,  extremely  important  to  insure  the  desired 
dual-frequency  behavior. 
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Analysis  of  a  circular-disc  printed-circuit  antenna 
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Abstract 

The  input  impedame  ami  the  electromagnetic  fields  of  a  cncuUr-disi  printed-circuit  antenna  are  found.  The 
antenna  is  driven  at  an  aihritiary  point.  Hie  present  analysis  lakes  into  account  t lie  radiation  through  the  side  wall 
ami  the  lunge  held  at  the  edge  ol  the  printed-circuit  antenna. 


I  Introduction 

Prnited  ciicuil  antennas  operating  at  microwave  frequencies 
.in  n>*.v  comnioiilv  used  on  many  (lying  objects  ami  satellites.  I'lg  I 
shoes.  i  smide  tin  nl.ndiM  antenna  etched  on  a  inn  rowave  punted- 
i  la  mi  hoard.  (lie  antenna  has  a  radius  equal  to  a.  and  is  driven  l>\  a 
•  oasi.rl  line  (Inough  lh<‘  gioiind  plane  at  ladial  distance  /*  Imm  the 

I I  nil''  I  he  thickness  and  the  dielectric  constant  ot  tie*  suhsliale 
o|  the  punted  t  to  ml  hnaid  .lie  denoted  </ and  respeelively  I  he 
peimeahilily  p  ol  the  siihstiale  is  equal  to  that  ot  ihe.m,  i.e  /;  p,, 

In  ptevuuis  stmlies,1  1  analylic.il  loinuilas  loi  the  elect mmagiiefic 
In-Ids  inside  the  aiilctuta  are  derived  hy  treat  my.  tlw  antenna  as  a 
resonant  cavitv  hounded  above  ami  below  by  conducting  plates  ami 
on  (lie  side  by  a  polled  magnetic  wall.  It  has  been  shown  that,  from 
this  simple  model  o)  the  antenna,  (he  radiated  fields  can  lie  obtained 
m  simple  analytic  bums  which  agiee  quite  well  with  the  measured 
data.  However,  this  simple  model  will  lead  to  minute  input  impedance 
and  infinite  electromagnetic  Helds  at  the  icsonant  frequency/  This 
difficulty  is  avoided  in  previous  woik  by  using  equivalent  circuits. 1 
The  limitation  of  this  approach  is  that  the  dielectric  substrate  must  he 
very  thin. 


(  tr<  uhr due  pnntrJ'Vircuit  antenna  with  coaxial  feed  at  r  -  h. 

In  the  present  analysis,  the  antenna  is  modelled  as  a  cavity  with 
limie  admittance  wall  on  its  side  surface.  The  admittance  surface  is 
determined  hy  the  radiated  power  and  fringe  Held  on  the  edge  of  the 
•tr?ferui,i.  Using  this  model,  the  input  impedance  can  he  calculated 
directly  from  the  basic  definition. 


/*■//*<  /■  .V  »  >4  t.f  first  received  JSth  %Ujfth  unit  in  revised  fnrrti  list  August  l ° -  ^ 
V?r/>  a  -.vith  the  Department  nf  t  tre tricot  I  nqineerina,  Hnivcnity  of 
ff'wston,  Hnt/shm.  h  \<j\  *7004,  USA 


2  Fields  inside  model  cavity 

Assume  ili.it  A,./  I  wlivn*  A,  w*  being  the 

augulai  tiequencv.  v>  lb.il  Helds  do  not  v.ny  m  tin*  ;  duett  ion  f’be 
("tal  lieldsdin*  t<»  .1  .  oiisl.iut  current  /„  main  tamed  by  the  coaxial  line 
at  r  h  .nc  given  ,o  f  i »l l« ■  ws 
In  region  I . 


1 , 

^  .  1  „./„  (A  |  r)  cos  n1,'1 

(la) 

Hr 

^  o An  Jn  (A  ,  r)  sin  n<J/ 

copr  n  -<i 

(Ift) 

/A,  • 

N  In  A,  (A  \r)  cos //<,*» 
l.j/j  „*r., 

(If) 

and.  in  region  II, 

/■/, 

y  ( A !  r)  t  f  A , r )  |  cos  nip 

Ca) 

Itr 

A  )r|/f,iy„(A,/|  +  C„Y„(ktr)\  sin  no 

u>vr 

Oft‘ 

'  /  |/?„./„  (A  i  r)  +  Cn  Y„  {k , r)|  cos 

U)|U  (,  .1 

(» 

where  ./’,  </nd  aic  liessel  functions  of  the  lust  and  the 

second  kind,  .iml  their  deitv, i fives,  respectively.  .1*.  H„  and  C„  a»e 
iousianis  to  he  determined  latei.  The  time  factor  <*/wf  is  assumed 
thrmighoin  the  an.iK  ms. 

The  following  houmlary  conditions  arc  then  imposed 

H.,„ 

r,„/,„  (air  al 

(da) 

ft,.,  in  It, i  (at  r  ft) 

(dft) 

(/,.,/>) ft «*.•>»  (atr  -  ft) 

(d<) 

Note  that 
that  the  c 

the  boundary  condition,  cqn.  }a.  is  equivalent  to  stating 
aviiy  is  bounded  atr  a  hy  an  admittance  surface  ytn.  The 

value  of  will  be  determined  latei. 

The  constants  An,  H„  and  C’„  in  cqns.  I  and  2  are  solved  to  satisfy 
cqn.  3,  with  the  following  result- 


A„  ~  ,!  {-/«r*i*)tn:r*i«)+/riA- >'«(*, «)l 

Z(1  +  h„)  f)„ 


-  />„  y„a,h)i 

/AitOft  I 

2i\+K)  n„ 

!  f/nhJfl 

:<V  os„) 


(4a) 

C  (A ,  />)  |  Y'A*  m)  +  n .  .»•„  (* .  a)|  (4ft) 

V„(A,ft)  (4c) 


whole  ftn  is  /oro  for  n  >  0  and  is  ci|Ual  In  I  lor  n  =  0.  (i  =  {vie \  )‘ 
Also, 


!}n  C(Aia)  t  /f|  l',..C<*|0)  (A) 

Nolo  lluii ,  will'll  k ,  a  is  cipial  to  otto  of  the  /cuts  o i  J„  (,v)  am!  wliott 
i„,  is  m-flocled,  /)„  will  vanish  am!  A„  and  ll„  will  become  infinite. 


3  Radiated  fields 


I  lie  radiated  fields  may  be  calculated  from  the  aperture  field 
at  r  -  a.  for  the  index  n  (called  the  nth  mode)  the  following  formulas 
aic  derived  for  the  radiated  fields: 


t-nOi" 


sin  (A nJ  cos  (?) 


cos  0 


—J’n{k„tt  sin  0)  cos  (6a) 


,.r  ,  I  Sill  |A| ,‘J  ios  ; 

'  ./.II  /  ««/  ,  Mil  W)  SHI  «v 

r  ktt  mu  ! 


/  //„./„  (A  t</)  t  ( K„(A  |  )  ami  A„  * 

4  Surface  admittance 

’I  lie  surlacc  admittance  y„n  associated  with  the  nth  mode  is 
ilc lined  in  eipi.  .ki.  The  real  pail  of  g,,,  g,,,  *  /b,„  *s  relalcil  t«»  the 

total  radiated  power  VTait  as  follows. 


4  // 

where  the  integration  is  over  iliccyhmlne.il  surface  at  r  -  a  hounded 
hy  0  The  radiated  power  may  he  obtained  by  integrating  the 

I’oyntmg  vector  over  a  large  hemisphere  using  eqti.  6.  The  result  is  as 
follows 

v.„  [d  *■*„)(*.,«»>-'/,  +  »r(l  '>n)/J|  <*) 

-"U  ' 

«!»•'.'  U  ' 

/i  |  |7,,  U„j  mii  «)|'  sill  «</0  (')a) 

and 

*  cos*  W 

/:  |  sin  0)1*  <//J  (0/») 

•  <>  sill  0 

In  the.ihuvc  lommlas.  A„»/ has  l»ecn  taken  in  he  inuch  less  than  units 
s<»  i hat  in  etpi.  r»  sin  (A„</  eos  0 )  is  appioxunatcd  l»v  (A0</  cos  0). 

I  he  imaginary  pail  o!  r,„  is  m.inil v  due  to  ihc  hinge  held.  lo 
estimate  it,  let  e«|it.  S  he  wiitten  in  teimsol  real  and  imagmaiy  parts 

I’.,  I  /„  ( A.  |  <v )  ( i ,(A  ,(/)|  1  /ii  ,(A  |tf  I  (Hi) 

In  i  lie  /emth-onler  llicmy,  I  lie  nulunum  is  nojileeleil :  hence  n,n  (). 
Hi.'  1 1 injie  t iel.l  is  jlsu  nejilee toil  lesiiltinyr  in  />„„  0.  The  resonant 
liei|ucncv  is  then  (leteinnneU  sini|>ly  hy  t lie  the  vth  zero  of./,’.  In 
:i  leeent  papci,'  I  he  etleel  ot  the  Innye  I  iel.l  is  taken  into  account.  It 
o  shown  that  the  el leet  of  the  Irinye  llel.l  is  to  lower  the  resonant 
t ie.|ucnCy  hy  a  factor  ei|iial  to  ( I  f  A)  1  ’ ,  where 

A  -  ~'k"  I  In  j"l't  f  I -772(1 
ire,  a  I  \ 2</ 1 

I  hits  the  model  cavity  will  he  resonant  when 

A,j  /,'„/(  I  l  A)1 til) 

l  snip  this  icsiilt.  it  is  clear  lh.it  /)„  should  vanish  when  e,„  0  and 

A  |ii  satislies  ei|n.  (  I  In  other  words. 

><„,  (».>/(?.  /., ( .» ) |  (12) 

" nh  «•  I  l  A)1  ’ . 

Ihis  completes  the  derivation  ol  the  loimuias  tor  the  internal  and  the 
ladiafed  tields. 

5  Internal  electromagnetic  fields 

In  practice  .1  ciieular  disc  prmtcdcircuit  antenna  is  opeialmg 
neai  I  lie  icsonanee  ol  the  n  I  mode  since  it  yields  a  strong  radiation 
m  the  j-ditecljon.  In  the  design  ol  such  an  antenna,  electromagnetic 
tiHds  aie  assumed  l<*  he  /ero  lor  all  oilier  inodes.  The  picscnt  llieoiy 
enables  us  to  check  the  Validity  of  this  assumption,  l  or  example, 
consider  a  prinlcd-cucmt  antenna  ol  modeiafc  thickness  with  J 
lH(»cm.  a  I  HXcm.  f  1  2-47«„  and  /»  (l*Mcm.  I  tom  eqn.  II. 

the  tesoiiaut  (requeue)  tot  n  I  mode  is  2-XVK»H/.  fig.  'til  shows 
.1  plot  ol  ./•,!  as  a  hmcliori  ol  r  lor  the  Inst  I  !  modes  in  the  antenna 
at  this  ic  sou  ant  t»cquenc> .  I  lie  fields  aie  obtained  Irom  eqns.  I.  7 
and  4.  v\  1 1 f ■  r,„  set  equal  to  /ero  lor  all  n  except  n  I.  The  value 
is  y.aku!a'cd  Irom  e«ins.  K  and  I  7.  h  is  noted  that  m  I  te.  Jli  t  the  li.dd 
it  the  edec  lor  the  //  I  mode  is  at  least  dH  lurhei  than  that  «t 
my  othci  mode  lhu\,  11  is  a  good  approximation  to  neglect  r„„  i«»r 
ill  other  modes  -  scept  lliat  o|  the  n  I  mode.  figs,  .'(u)  and  '(m) 
d'ow  the  fields  m  tin-  antenna  at  frequencies  considerable  lowei  than 
die  •eson.mi  fro* 1  w‘\  The  n  I  mode  is  'till  dominant  in  the 
.a\ils  jt  these  tie. men. t-  ' 

5  Input  impedance 

I  h.*  input  unpt'damc  udehm'd  as*' 
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where  h\.  and  Wm  arc  time-average  stored  electric  and  magnetic 
energies. 


and  Pr  -  I‘raa  *■  I'/,  +  with  Pa  ami  /'i.  being  the  power  lost  in 
dielectric  substrate,  and  that  lost  as  ohmic  heat  on  the  conducting 
plates,  respectively. 

Assuming  that  the  n  ■■■  I  mode  is  dominant,  the  above  quantities 
can  be  expressed  in  analytic  lonns: 

1' 

Pr, i,r  ;  rri/4 / 1 H (  J ( ( k  t a )  t  (  ,  1 1  ( A I  fl  H  e,i 

.\.)h 

2t><  lan  A  It', 

k  j»J)> i!!'" 

HA  , 

i  >/t ,  r  |  v  *’  t ./ ,  j„  j.  » | 

*  f  n»:s 

i  t(H,c;  t  «;c,)|e'i’y,  k,  ./„y.  JiK,, |j;s } 

.  .  5 1  el  1  1 2  I  V  *  ( y  o  -J'l  |ii|h 

I 

+  +  J\ ) 

*  *■  y ;> 

+  1 y,y, 

k/.k,  J  y,  k,  •  l7,K,  |J;Si 

where  |/  |»  means  /M)  Jili).  it,  ~  (rrqri/  )  1  ;  is  the  skin  depth  ot 
flu-  ci inducting  plate  at  the  operating  frequency,  o  is  the  conductivity 
ot  the  plate,  and  Ian  ^  is  the  loss  tangent  of  the  diclcctiic  substrate. 

Hie  above  formulas  are  deliver!  tiom  indefinite  integrals  involving 
Mcssol  functions.  Tliese  formulas  are  listed  in  Appendix 

I  v  l  shows  theoretical  input  impedances  of  circular  disc  printed 
*  m  ml  aiPcnn.n  with  dlJteicnt  driving  points.  Agreement  with 
measur'd  data  is  seen  to  he  good.  Note  that  the  antenna  was  etched 
■  m  a  rather  thick  printed-circml  hoard  (>//</  OOMS)  (he  input 
resisiaui •'  al  resonance  agrees  with  the  measured  value  beltei  than  the 
previous  theory  *  docs. 
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9  Appendix 

In  deriving  lire  formulas  lor  It',,  and  irm,  the  lollowing 
mieg'.il  birmiiias  are  used : 

1  '[y.nti'./v  S  :jy..u)|;  *■  |y,i v»rJ  1  <i4i 

•  V 

I  v  | ./ ,  (v)| :  ilx  -  (  :<  V-  *  l)|./„(<r)|-’  +(.v  I »(./,  (,v  >|  -  } 

(I'M 

iviy.i.viivjt  =  \v'<tj.txtV  ./..iviy.io:  not 


I  1  /,  1  v)  r,t  v)./v  j  |y„(.v) » „(v)  +  y,  i.v)  k,  oh  ( i«» 

•  V 

|  v./;(v|V,'(v)Jv  -  (I  +  l.r:)((y„(.rl>'„(jr)  +  y,0)>i(v)| 

4  i.v:py.(.viy.|.r)  y,(Af))-,(.vi 

y,(.v))'i  u)| 

In  n|iis.  M  lf>7  may  be  replaced  by  V.  These  formulas  arc  derived  or 
directly  quoted  from  those  in  Reference  7,  pp.  1.14  !  Ms,  and  trom 
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Radiation  patterns  of  a  circular  rnitrosirip  antenna  have  been  calculated  with  Green’s  functions 
for  a  two-layer  medium  and  compared  with  previous  results  which  are  valid  for  thin  slabs  and 
small  dielectric  constants.  A  mode  expansion  technique  is  then  used  to  find  fields  and  resonant 
Frequencies.  The  theoretical  results  for  the  input  impedances  are  found  and  are  in  good  agreement 
with  experimental  data. 
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I.  INTRODUCTION 

A  microstrip  antenna  consists  of  a  metallic  patch  print¬ 
ed  on  lop  of  a  dielectric  substrate  over  a  highly  conducting 
ground  plane.  To  calculate  the  radiation  characteristics  of 
such  a  structure,  we  must  use  the  field  formulation  fora  two- 
layer  medium,  which  has  been  a  problem  of  historical  inter¬ 
est. '  '  For  a  circular-disk  microstrip  antenna.,  theories  have 
been  developed  by  using  frec-space  Green’s  functions. In 
this  paper,  we  first  calculate  the  radiation  patten  by  using  the 
two-layer  stratified  medium  formulation  and  compare  with 
previous  results'1  which  are  known  to  be  valid  only  when  the 
dielectric  constant  of  the  substrate  is  not  high  and  when  the 
substrate  is  thin  as  compared  with  the  radius  of  the  disk.  Wc 
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next  compute  the  resonant  fields  by  using  a  mode  expansion 
technique.'  K  We  show  that  the  results  obtained  for  the  input 
impedances  agree  very  well  with  experimental  data.1' 


II.  FORMULATION 


Consider  a  microstrip  antenna  modelled  as  a  current 
sheet  placed  at  the  interface  of  a  dielectric  slab  on  top  of  a 
perfectly  conducting  ground  plane  (Fig.  1 ).  The  current 
sheet  can  be  represented  by  a  current  density 

J(r)  =  K(x,  yy>{z).  (1) 


The  electromagnetic  field  components  due  to  a  point  source 
radiating  in  the  presence  of  a  stratified  medium  have  been 
formulated  in  integral  forms.'  By  decomposing  the  current 
shed  into  its  jc  and y  components,  we  can  show  by  superposi¬ 
tion  principle  that  the  field  due  to  the  source  given  by  Eq.  ( 1 ), 
with  lime  dependence  e"“  omitted,  is 


£,  =  f  d  V  [cos(<#  '  +  X  )K*  +  sin(<*  '+*)*,] 

j  OlT(t)€  J 

<x  J '  dkr  k2„(\-R 1  V  ’(*„  I  p'  -  p  I ),  (2) 

H,  =!~  f  d  V[sinW  +  *)*,  -  cos(<* '  +  x )K, ) 

j  8  7T  J 

;X  [  dkr  ^  ( 1  +  R  TV  7/  V  ’(*„  |  p'  -  P  | ).  (3) 


wlierejk,  —  (k  1  -  k  /?1M  and  R  "are,  in  this  case,  the 
reflection  coefficients  of  the  two-layer  stratified  medium  and 
aregivenby  (A3)  and  (A4),  and //'p''(Jc(,  |p’  -  p|)  is  the  first- 
kind  Hankel  function  of  the  first  order.  We  note  that  the 
dependence  of  the  field  due  to  Kr  is  just  that  due  to  K,  rotat¬ 
ed  by  90°. 

To  compute  the  radiation  pattern,  we  arc  only  interest¬ 
ed  in  the  limit  when  |p  p’|  -*«  and  z-  -»oo.  In  such  a  case, 
tfi '  +  V  =  and  Eqs.  (2)  and  (3)  become 


/:.  ~  -  —  I  d  r  (cos^A',  -f  sin^A.',) 

8  note  J 

X  J’  dkf,  k2\\  -  R  ,My‘ ip’  -  Pi). 


(4) 


//,  -» — -—  fr/Y<xin^.  —  cos 4>KS) 
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X  f '  dk,  -~(i  +  R  "V*  'll \'U„  Ip'  -  pi).  K-  EJ :  (—  p)  cos n*. 

J  .  k.  t  / •) /u\  \  a  / 


The  other  components  of  K  and  1 1  field'  can  he  derived  from 
Eqs  (2)  and  (3).  I)y  using  the  steepest  descent  method,  it  can 
be  shown  that  the  electric  field  components  in  the  far  field 
are  given  by 

/•;.  - —  k  !  sin 0  cos 0  1 1  -  R  ,M(0  )|  — 

4  mat  r 

X  J  <f'r’(cos<M\  +  sin^A^K-  (6) 

Er  '  k  7  cos  0  1 1  -  R  ,M(0)|  — 

4rr<t)e  r 


J  </ ’r’ [cos^Af, 


+  sin^A',  )e 


E,  --=  ^-[l+*"(0)]£- 
4rn  r 

X  J  cf  V(sin^/f„  —  cos<fiKy)e  .  (8) 

In  the  above,  we  have  omitted  the  singularity  contributions 
in  the  steepest  descent  approach,,  which  can  be  shown  to  he 
negligible  (see  Appendix  A).  We  also  notice  that  the  far  field 
approximation  of  the  field  is  similar  to  that  due  to  a  current 
sheet  in  free  space  except  now  that  the  TM  waves  are  scaled 
by  a  factor  1 1  R  IM(0)]  and  the  TE  waves  by  a  factor 
1 1  -4  R  "(0) ). 

III.  RADIATION  PATTERNS  DUE  TO  A  CIRCULAR  DISK 
ANTENNA 

To  find  the  radiation  pattern  due  to  modes  excited  on  a 
circular  disk  antenna  printed  on  a  microstrip  board,  we  need 
to  know  the  current  distribution  on  the  disk  for  various 
modes  The  approximation  to  the  current  distribution  can  be 
obtained  by  studying  the  modes  of  a  circular  cavity  sur- 
tounded  by  a  cylindrical  magnetic  wall.  For  such  a  cavity,, 
the  field  inside  is  given  by"’ 

E,  =  E„  J„  [( ft„m  /a)p  J  cosn.0,  (9a) 

,,  donm(  |/r  //)„„,  \  ,  /nu\ 

H>,  =  —  ~y-  —  E«  Ji  I - P ) Sinn*’  (9b) 

(P„,„/aYp  \  a  ) 

H.  -  E"  J-(  —p)  cosn*’  <9c> 

where 

.  ft, . /|rr(//,r,)'/?|. 

Since  the  tangential  II  field  has  to  vanish  on  the  mag¬ 
netic  wall.  //#(r  0)  -  Oory  '(//„,„)  0,  where  a  is  a  radi¬ 

us  of  the  cavity.  The  current  distribution  on  the  cavity  wall  is 
given  by 

K  -  H,.i  +  Hf  p.  (10) 

1  hcrefore  the  current  distribution  for  various  modes  is  given 
by 

..  . f,»  ,.  ,  (P™  \  ,  ...  . 

kt  -  -  - - ,-E„J.\ - p  Isinn^,  (I  la) 

(ft.,.Ja)p  \  a  '/ 
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K  - — p  cos«<».  (I  lb) 

(0,,,.,/tf)  V  a  / 

We  cun  assume  F.q  (I  I)  to  be  the  zeroth  order  approxima¬ 
tion  to  the  current  distribution  of  a  circular  microstrip  an¬ 
tenna  and  use  it  to  calculate  the  radiation  pattern. 

Substituting  Eq  (I  I)  into  Eqs.  (7)  and  (8)  (see  Appen¬ 
dix  H),  we  obtain 

t$k  r 

E„  -  <.<//«  .//i,)"-’^— (  -  ')"  ’  1 

2  r 

(ft . /a) 

/  /■-„  cos 0  cos n<j> - - ; - - — — 

( /l, ,,,,/ay  —  k  ■  surf) 

A  I aJ„  (  ft.,,.,  )J (ka  sniff )  ]  [  1  -  R  1  M(0 ) ) .  (12) 


..  ,  ,  E„  sinnd) 

*  2 r  (ft /a)k  sinf) 

y.JAft . )J„(ka  sin0)l  1  +  R  "(0)].  (13) 

When  n  -  I,  the  nuxle  has  the  lowest  resonant  frequency 
and  has  been  extensively  studied.  In  Figs.  2  and  3,  we  plot  the 
radiation  field  when  to  —  ut„,„  for  different  e,  =  (f  ,/e)  and 
d  /a  ratios.  The  electric  fields  are  normalized  with  respect  to 
rE,pJ„(  ft.,,.,)-  Long  cl  at.1'  calculated  the  radiation  field  by 
assuming  magnetic  current  sheet  on  the  side  wall  of  the  cav¬ 
ity  and  using  free-space  Green’s  function.  Their  results  are 
plotted  in  dotted  curves  for  comparison.  We  note  that  when 
d  /a  is  small,  both  result  agree  very  well  except  for  0~9O". 
Since  their  formula  does  not  account  for  the  dielectric  effect 
of  the  substrate,  we  conclude  that  the  dielectric  effect  on  the 
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radial  ion  field  is  small  when  d  /a  is  small  However,  when 
i /  /,i  is  large,  noted  discrepancies  are  observed.  These  dis¬ 
crepancies  aie  larger  when  t,  is  increased.  We  find  that  the 
radiation  field  increases  in  intensity  when  d /a  increases. 
Also  the  directivity  in  the  normal  direction  improves  when 
d  A;  increases  and  f,  decreases  Whetu,  is  large  or  d  /a  small, 
more  power  is  directed  near  the  horizontal  direction. 

An  additional  clfect  of  the  dielectric  is  the  presence  of  a 
guided  surface  wave  I  hat  carries  power  only  in  the  horizon¬ 
tal  direction  al  0  W.  However,  when  the  microstrip  dtsk 
is  mounted  on  a  finite-si/e  ground  plane,  this  guided  wave 
can  be  diffracted  into  other  angles  and  gives  rise  to  radiated 
power  When  the  dielectric  layer  or  the  ground  plane  is  loss- 
1  v,  the  guided  wave  has  a  complex  wave  number  [see  Eq. 
(Afi)l  In  such  cases,  the  guided  wave  decays  exponentially 
away  from  the  source  as  its  power  is  dissipated  in  the  dielec¬ 
tric  layer  or  ground  plane 


IV.  MODE  EXCITATION  OF  A  MICROSTRIP  DISK  BY  A 
PROBE 

In  the  limit  of  small  d,  the  microslrip  disk  acts  like  a 
high  Q  resonator  and  the  field  inside  the  disk  resonator  re¬ 
sembles  that  of  a  closed  resonator  with  a  magnetic  wall  on 
the  side. The  modes  of  the  resonator  cavity  can  be  divided 
into  the  natural  modes  and  the  unnatural  modes, ’  "  in  terms 
of  w  inch  arbitrary  K  and  II  field  can  be  expanded.  We  write 


K  1  <’ . rK . r  f  X  . - 


<  1 4a) 


II  ^  "f  ^ 

tt.m,  p  fi,m,  p 

where  n  denotes  t  he  <t>  variation  m  the p  variation,  and p  the  z 
variation  of  the  field.  In  the  above,  j  |  forms  the 

natural  mode  and  V„mp  is  the  unnatural  mode  for  the  electric 
field,  and  G„,„,,  is  the  unnatural  mode  for,  the  magnetic  field. 
The  natural  mode  fields  satisfy  the  following  equations: 


^  X  E*m/,  —  knmpHnmpt  V  X  —  knmpEn  mp 


<v‘  +  *:,Xr)=° 

and  the  boundary  conditions 


(15a) 

(15b) 


(15c) 


n X  E„mp  =  0,  on  S,,  n-E„mp  =  0,  on  S„ 

«  H„1(,=0,  on  Sr,  n  X  =  0,  on  S„, 

where  S,  indicates  the  electric  conductor  surface  and  S„  the 
side  magnetic  wall.  The  unnatural  mode  for  the  electric  field 
satisfies 

(16a) 


k  K 

nmp  nmp 


Vt K 


where 

(V2  +  k  L,r)4’n,„P  =  o 

and  the  boundary  conditions 
'I'nn.r  =  const  I,  on  S, 
=  const  2,  on  S, 


(16b) 


and 


dp 


=  0  on  S,, 


(16c) 


where  Su  is  the  top  conducting  disk  of  the  cavity  and  S,  the 
bottom  conducting  disk.  Since  V-F„m/>  =  (1  /k„„,p)V2if>rmp 
=  -  we  see  that  if  k„mp  *0,  the  electric  field 

corresponds  to  the  field  due  to  static  electric  charge 
distribution^,  =  -  € , k, „„r\l>nmp  permeating  the  cavity. 

The  unnatural  mode  for  thejx>agnetic  field  satisfies 

k„mpk* . .  =  V  (17a) 

where 

(V-  4  k  w . ,-0  (17b) 

and  the  boundary  conditions 
<t> . .  0.  on  S„ 


and 


=  0, 

dl 


on  S,.. 


(17c) 


We  find  that  V-G . .  =  —  k„n,r<f>„mp,  which  implies  that  the 

magnetic  field  G,„„(,  corresponds  to  field  due  to  static  mag¬ 
netic  charge  distribution/^  =  -  pik„„,n<t>„,„n.  Since 


Vx 


the  unnatural  modes  are  also  termed  the  irrotational  modes. 
To  solve  for  c  ,  h„,„  ,  we  make  use  of  the  identity7 


1 


V-(Ex  H„,„)  dV 


i 


(VxE)-H„mpdy- 


l 


(VXH)  EdV 
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n  X  EMI . ,,  dS, 


(18a) 


(lie  relations  V  >  -  k . ami  V  >  K  im/t  ,11, 

anil  t  lie  mode  orthogonality  properties  to  obtain 


J  ’<  '  EMI... 


dS. 


Ui/fl  \hnmp  kttlll/1l 
Similarly,  using 
V  (K . .  II )dV 

j  (Vx  F. . „•!!  -  V/H  dV 

n  x  H  E ( IS 

mil  the  Maxwell  equations  V  x  H  -  iiot ,E  4-  J, 


1 


(18b) 


(1 9a) 


v  >  E . . 

=  H . .  we  arrive  at 

^  i  ••»/>  ^  »»»»(/• 

->  <>-*.<• . .  -  J  J-k . .  dy 

|  >'  ■  H-E . .  dS. 

(19b) 

1  o  solve 

for w  e  make  use  of 

1 m  7 

» .  »• . „.-V  /U)dY 

J  »  X  H-E . .  dS 

(20a) 

to  obtain 

‘c>(  1  J . ,, 

J  J  E . .  di'  J  n  X  H-F„,„(,  dS 

(20b) 

for  X . . 

.  we  make  use  of 

J> . 

VxE-  E-VxCJ . r)dV^  |  n x E-G,„ 

„rdS 

(21a) 

to  get 

. .  »XK •(.' . .  dS. 

Js 

(21b) 

I  quations  (1Kb),  (1%).  (2()b),  anil  (21b)  can  be  solved  ap¬ 
proximately,  if  hi  is  such  that  nr’/r,*  ,~A  In  such  a  case, 

we  get  K K„„,(„  II  l( . ,, .  Hence  the  left-hand 

side  of  ( 18b)- (2 lb) can  be  approximated.  Wc  have  from 
( 1 8b)  that 


ri  x  EMI . .  dS~Z,h . .  | | '  dS,  (22a) 

where  is  the  surface  impedance  of  the  metallic  surface 
given  by  Z  (1  i)/o< 5, ).  Since  the  real  power  dissipated 

on  the  metallic  surface  is  given  by 


and  the  Q  due  to  this  loss  is 
toft ,  J 


Q,„ 

we  have 


n„,„r\-dy 


2  P, 


<•>!<  i 

2  P,, . 


ii  X  E-H,,,,,,,  dS*  h . „(  1  -  /).  (22b) 

Js  V,/ . 

From  (19b),  we  can  simplify  the  right-hand  side  as 

_[  "  *  H  E . .  dS~  Y,„e . ,,  |  |  E„,„p  |-  dS,  (23a) 

where  )’„,  is  the  surface  admittance  of  the  surface  -S’,,, .  In 
essence,  we  have  replaced  the  surfaces1,,,  with  a  lossy  mag¬ 
netic  wall.  In  general,  Y,„  0',„(1  in),  where  «  is  nonun- 

ily.  Similarly,  we  can  relate  (23a)  to  the  Q due  to  the  magnet¬ 
ic  wall  loss  and  have 


J  "  H  E . .  dS-  . ,,(l  -  ia),  (23b) 

where 


& . 


!»•: . r[dV  ~  ruf, 

2 P,  ~  TP~ 


(23c) 


is  the  Q  associated  with  the  radiated  power  loss  P,  of  the 
nmp  mode. 

Therefore  ( 1 8b)  and  ( 1 9b)  can  be  solved  for  c„„,,,and 
h . .  giving 


'cy/^l  4  (I  fi)/()„ . j.l,  dV 

k  i.„:  -  ">:p  if  I  [  1  +  ( 1  *  ]  1 1  +  («  +  i)/Q, . .  ] 

^  ^ I  r  dY 

' . . ”  ’  «’  ><7.1 1 1  <  1 7'»w~]  [iT(7Tm2~T 

Noting  that  the  second  term  on  the  right-hand  side  of  (20b)  is  small,  we  have 

-  - 1  (  J  i . .  dy 

Hoe |  Ji 
and  also  t hat 


(24) 

(25) 


(26) 


8 *0.  (27) 

Wc  shall  assume  that  the  current  source  due  to  the  probe,  when  d-^A,  is  given  by 

J  :IS(x  b)6(y).  (28) 

W'c  sec  that  due  to  the  uniformity  of  the  current  on  the  probe  along  z,  only  the  E II . and  F,„„„  modes  will  be  excited,  that 


I'llO 
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J  z/Mx  -  l>  y>(  y).  (28) 

We  sec  that  due  (ol  lie  iinifuimily  oil  he  current  on  the  probe  a  long  c,  only  the  K . II . and  l’„ . modes  will  be  excited,  that 

is,  only  modal  tields  with  no  z  variation  are  excited.  Since  the  current  probe  has  uniform  current  distribution,  there  is  no 
s  ha  rye  aecuimilation  in  the  cavity.  Thus 

if  (..in  be  shown  that  no  F . .  modes  are  excited  Similarly,  since  all  the  TIi  modes  of  the  cavity  have  z  variation,  none  of  them 

aie  e  veiled  As  such,  the  natural  mode  of  the  magnetic  wall  cavity  which  is  of  interest  is  the  TM  mode,  whose  held  is  given  by 

K zE . V„  | ( /i /a)fi  j  cos«<4,  (29a) 

vvhcie  J  ‘  (  // )  0  We  have  assumed  that  the  probe  is  placed  at  <#  0  so  that  K,„„„  is  an  even  function  of  d>.  The  magnetic 

held  of  the  natural  mode  can  be  derived  using  ( 1 5a)  giving 

/  u  \  l 

(29b) 


(30) 


- J„  pj  stn/id  4-  6J  ’  ^  — —  f> j  cos n<t> 

lali/.ing  constan 

I  |K, . I’ r/r-  111 . r’ dr,,.  1. 


(P . /a)f. 

In  the  above,  /:,„„„  is  a  normalizing  constant  such  that 


Thus  E 

E, .  -  (2/nd) 


can  be  derived  easily. 

(IJA,,.,)1^ 

ujjfi . )  (PL 

Also  A: . in  Kq  (15)  is  given  by  (P,„„/u)  in  this  ease.  Using  the  above,  we  find  that 


j  ,)•!■: . dr  /■; . j,,^1 .  i> )/</ 


I  Inis  we  conclude  that 

""/',/d  1 1  I  (I  I  .  . J.,\(P, ...,/“)!> 


I  :~i  '  -  \  . .  ---  J 

(P . /ay  w/i ,f , [ I  I  (I  I  i)/Q., . J[1  +  (</  t  >)/Q,  | 


'-(~f 


cos  m}> 


(31) 


(32) 


(33) 


and 


II?- 


Idk . a;;! . J„[(P . /a)h\  _ 

(  P„„,/a)-  -'«>>,  p  4(1  4  /)/(?„ . ][l  t  (a  f  i)/Q, . ] 

p - " - -  siting  f  4>J '  p )  cos n<t> 

.  (P . /a)f,  V  a  /  \  a  / 


(34) 


when  k  . ,e  ,.This  is  true  when  Q,,  >1  and  Q,  >  1 . 1  he  resonant  frequency  of  the  disk  is  given  when  the  denomina¬ 

tors  of  (33)  and  (34)  vanish.  Thus  the  resonant  frequency  for  the  nnw  mode  is  given  by 


_  ?■-  If,  +  +2_±_£')|  ’  ■ 

II  a...  A  e.  _  !\ 

(\  *_LL  " 

)r:'«  V  ~  2 Q, '  2Q. ) 


(35) 


-  P.z,  /,  lj_/_  rr 

I  a  V  2Q, .  2Q, 

I  he  frequency  shift  from  the  magnetic  wall  model  is  due  to  real  and  reactive  power  absorbed  by  the  metallic  wall  and  real  and 
reactive  power  leaking  out  of  the  disk  resonator.  Qt,  can  be  computed  since  we  can  compute  P<t  ,  the  power  dissipated  on 
the  metallic  surface  due  to  the  in  no  mode 

We  can  also  compute  (2,  through/',  .  the  real  power  radiated  by  the  disk.  Wenotclhat(l  -  nr}P,  is  the  complex 
power  leaking  our  of  the  disk  resonator  through  the  side  apertme.  This  power  includes  the  power  in  the  guided  surface  wave 
and  spherical  wave  radiating  into  infinity.  Thus  this  (rower  is  also  dissipated  in  the  dielectric  loss  and  copper  toss  outside  the 
cavity,  Che  reason  there  is  complex  power  leaving  the  cavity  i>  because  of  the  actual  nonvanishing  tangential  H  field  at  the 
aperture  of  the  cavitv  To  find  a  first -order  approximation  to  this  nonvantshing  tangential  H  field,  we  assume  a  current 
distribution  on  the  disk  due  to  a  natural  mode  given  by 

(36a) 
(36b) 

l ’ving  the  st  ratified  medium  formalism,  it  can  be  shown  that  this  current  source  gives  rise  to  nonvanishing  tangential  II  field  at 
the  aperture  given  by 


E, 

n 

p  )  snind. 

""/'.  /' 
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-  j  „ 
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«i"  -w(^r  + * 

\Pip  Jn  kt,  smkud 


-**r 

Jo 


*„(*,.) 


cosk ,  ,  (2  +  d) 
coskud 


(I  +  R'")J’m(kl.p)dk 


) 


where 


*.(*„>  =  - 
<*.,<*„)=  - 


-~^~JA0.Mkra), 

Zk.wft, 

F  k  O1 

»»»n  »  /  />  .  >(/i  v 

—  )•/;(*,.  a) 


«W«l( /a)2  -*/.] 

Thus  an  approximation  to  P, _ (l  —  ia)  is,  from  (23b)  and  (23c),  given  by 


(37) 

(37a) 

(37b) 


P,  ( 1  -  tor)  =  -  i-  £  n  X  HE„ 


dS 


-  iJT 


ad 4>  dz  = 


. ft  L  } 


JAft.,m)FUo)0+S.m), 


where  F(co)  is  given  by 

h*i-  -*£-  f%^<> +*"><«„+  r», 

1  ft  L  2o  Jo 

Consequently,  we  find,  with  the  use  of  Eq.  (31),  that 

PI 


4i(op , 

[J  ;(*,,</)  ]:  tan£(i</ 


l (ft^/o)1  -  k],\  klt 


(1  +  /?1M)  dkp. 


a  +i  (PLY 

Q.  \ak,J 


rL 


F(to), 


(38) 


(38a) 


(38b) 


d  [  ft  L  -  "  . 

which  can  be  used  in  Eqs.  (33)  and  (34)  to  find  K  and  H.  The  integral  in  (38a)  can  be  integrated  along  the  Sommerfcld 
integration  path  numerically.  The  above  approach  is  similar  to  the  perturbation  approach  used  in  finding  the  resonant 
frequency  shift  of  a  microstrip  disk  resonator."  12 

V.  INPUT  IMPEDANCE 

After  having  studied  the  modal  excitation  of  the  microstrip  cavity,  we  can  compute  the  input  impedance  of  the  antenna 
quite  easily.  First  we  find  the  complex  power  P  delivered  to  the  microstrip  cavity  which  is  given  by 


--L 


EJ  *dV. 


Since  P  is  also  given  by 

r  =  z„\i\\ 

making  use  of  Eqs.  (28)  and  (33),  we  find  that 

Z  M11+0+  0/ga.„„  ]F-  LJ  l[{ft.m  /a)b  ) 

"  “  ( -ruW,[l  4(1  4  i)/Qd__  )  [  1  +  (a  +  i)/Q,_  ] 

The  above  can  be  simplified  to 

icopdJ  l[(ft„m/a)b  ]{P^/af[\  4(1+  i)/fi,_  ](2  -  <5,,„) 


Z„  =  - 


tr!(A,/a)J-«*P.€,l  1  +d  +0/Crf„„][l  +(<*  +  ‘VQr.m.)\(PL  -*Vi(A,) 


(39) 

(40) 

(41a) 

(41b) 


In  Figs.  4(a),  4(b),  and  4(c),  wc  plot  the  conjugate  of  the 
input  impedance  using  Eqs.  (41)  and  compare  the  results 
w  ith  experimental  data.'*  (The  conjugate  of  impedance  was 
plotted  because  we  have  used  e  in  our  analysis.)  We  find 
that  we  can  match  the  input  resistance  very  well.  Also,  we 
have  slightly  altered  f,  from  the  manufacturer's  specifica¬ 
tions  This  is  done  in  order  to  line  up  the  resonant  frequency 
with  that  found  by  experiment.  Since  f,  was  specified  to  be 
2.47  t  0.04  m  the  range  of  8- 1 2.5  GHz,  it  was  uncertain  that 
t .  remained  within  the  range  at  3  GHz.  Also  we  find  that  our 
approximation  to  the  reactance  is  not  as  good  as  predicted  by 
e  xperiment.  This  is  because  ,Y  arises  out  of  the  energy  stored 


in  the  near  field  of  the  driving  source.  Since  we  have  used  a 
single-mode  approximation,  this  near  field,  which  is  singu¬ 
lar,  is  crudely  approximated. 

From  Eq.  (41)  wc  notice  that  when  the  Q  of  the  resona¬ 
tor  is  low,  the  bandwidth  of  the  microstrip  antenna  is  larger. 
It  thus  seems  that  a  way  to  increase  the  bandwidth  of  the 
antenna  is  to  increase  the  radiation  loss  by  increasing  d  / a 
ratio. 
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I  I]  EXPERIMENT 


Fit;  4  Input  impedance  compared  with  experimental  data.  Theoretical 
curves  are  plotted  for  d  —  0. 16  cm.  a  -  lR8cm,f,  —  2.60,  n  -  1, 
a  -  5  8x10’  rnlio/m,  and  loss  tangent  =  1.8  X  10  (a)  b  =  1.88  cm,  (b) 
b  141  cm,  (c)  b  ~  0.94  cm. 


ami  Army  Research  Office  Grant  DAAG-29-78-G-0198 
with  the  University  of  Houston. 

APPENDIX  A:  INTEGRAL  APPROXIMATION  BY  THE 
STEEPEST  DESCENT  METHOD 

In  the  text  of  tins  paper,  we  can  simplify  integrals  (4) 
anil  (5)  provided  that  we  can  simplify  the  following  integrals: 

/,  I  «/*,.  yd  R  ,VV‘  7/;,"(A„||»'  pi) 

<AI) 


I  r  V-*»  !»"(*„  |p' -pi). 


(A2) 


The  inner  integrals  in  (4)  and  (5)  are  just  derivatives  of  /,  or 
/;.  For  a  stratified  medium  shown  in  Fig.  1, 


R 


IM  _ 


I  —  C 


2k ,  r/ 


I  f  R 


1  MJW|  <l 


2  sinA„t/ 

sin kt,d  4  H(,/t)(k,/k„)cmk,,d  ' 


(AJ) 
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I  U  i  A )  I  hr  siixpoM  desccni  palli  (SDP),  ihr  Matinnar)  point  ikrs ).  the 
hiiim  h  pomi  {k  ),  and  the  polc(AJV,)on  the  complm  kr  plane. 


»*■*""  -V,“ 7T7<* 

I  D  II 

I  -  K  e 

2  sinA.d 

= - - - ,  (A4) 

sinA„<f  +  /( p/p , )(A w/k, )  cask 

where  A1,.  --(A;  —  A  2 )1-2.  The  locations  of  the  poles  are  giv¬ 
en  by  the  transcendental  equations 


waves,  and  for  TM  waves,  there  is  at  most  one  pole.  The 
location  of  the  pole,  in  the  limit  that  d— »0,  is  given  by 

A,,.~A  +  l(A  V/VAJXA’  -  A;)2.  (A6) 

Physically,  this  implies  that  when  d  is  small,  there  is  at  most 
one  guided  mode  for  the  TM  waves  that  can  be  excited,  and 
for  IheTE  waves,  there  is  no  mode  excited.  This  is  true  of  the 
microstrip  substrate  since  the  dielectric  layer  is  thin. 

It  can  be  shown  easily  that  (A3)  and  (A4)  are  even  func¬ 
tions  of  A , . .  Therefore,  there  is  no  branch  point  at  A(,  =  A,. 
When  |p’  p|  »«>  we  can  replace  ’(A(>  |  p’  -  p|)  by 
(2/irA(,  |  p’  -  p|)l/V*,'"‘  M  in  (Al)  and  (A2).  Thus  there  is 
a  stationary  point  at  A(„  =  A  sin0,  where 

»  tan  '(|p’  -  p|/z). 

The  locations  of  the  singularities,  the  stationary  point,  and 
the  steepest  descent  path  passing  through  it  are  shown  in 
Fig.  5. 

The  pole  gives  rise  to  a  surface  wave  which  decays  rap¬ 
idly  front  the  surface  since  k,r  is  pure  imaginary.  Thus  when 
|p’-  p  |  *  oc  and  z  ♦  oo .  in  other  words.'when  we  are  in  the 
far  field,  the  dominant  contribution  to  (Al)  and  (A2)  comes 
from  the  stationary  point.  With  the  steepest-descent  analy¬ 
sis,  it  can  be  shown  that 


tanA  ud 


i:tiiA  ,,d -■  -  i 


for  TM  waves 


for  TE  waves. 


It  can  be  shown  that  the  above  can  only  have  solutions  in  a 
certain  domain  of  the  complex  plane,  and  thus  the  nature  of 
ihe  waves  due  (o  the  pole  contributions  can  be  predeter¬ 
mined.  Wlteiw/  is  small,  there  is  no  pole  contribution  for  TE 


1  ,,<*  a  i> 

- - 

i  I  r  —  p  | 

i  '  i* 

-  —  |1  -R  ,M(0)1  - - . 

/  r 

J  If  |» 

/.-  -  1 1  +«"«?)]  4 - Jr 

i  I r  —  P  I 

1  III* 

~  -  W  3  K  "((>)]- - . 

i  r 

We  can  use  the  above  to  derive  Eqs.  (6)-(8). 


APPENDIX  B:  EVALUATION  OF  THE  E-FIELD  FOR  A  CIRCULAR  DISK  ANTENNA 

First  we  shall  evaluate  the  A’*  component  given  by  Eq.  (8).  It  can  be  shown  easily  that  Eq.  (8)  can  be  reduced  to 

flkr  pt 

A.  -  -  1 1  +  R  M(f/ )]  - —  I  p  dp’  d<t> ’  [A*  cos(</i  ’  —  <b)+  K(,  sin($  ’  —  <t> )]  exp[  —  ikp'  sin//cos(<£  —  4>  ’ 

4  iri  r  J„  J „ 


Letting  the  integal  part  of  the  expression  above  be  /,,  substituting  with  Eqs.  (I  la)  and  (1  lb)  results  in 

/,  If  P  dp'  dd '  i((  ,/ft  - - -J,\———p'\ sin«d  ’  cos(^  ’  —  d )—  J (— '•—/>’)  coswd  ’  sin(^  ’ 

J"  j -  U0 /uV>’  \  a  /  \  a  / 

•  exp|  ikp'  s\nO cos(d  -  d’)|- 
Using  the  lad  that 

/(O  ^  j  ihl>  cvp<  m  cosei  inji  t  \inrr), 

andthat./  „(t)  (  I  )V„(v).  it  call  be  shown  that 

|2ir 

coswtfr  ’  sm(<£  ’  d  )  exp  |  ikp  sin//cos((A’  <k)\d<!>‘  «(  — J„{kp  silt//)  sin/n/i 

kp  silt// 


mtnk  ’  ci  >'te< '  if’  )  exp|  (/./■  tricos  (/’ '  d)l  dj>‘  2,r(  /)"  1  'J'„(kp  sin// )  siting. 


I  hciclori  'M2)  becomes 
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/,  -  i((  ,/p  E„n(  -  i)“  '  '2nsnm<t> 

Recognizing  the  integrand  to  he  exact  derivative,  we  have 

/, 


,  fu  V<"  | J  ‘  (kp  sinty )  J‘„  [( P„,„ /u)p‘  J J„ (kp'  sintf ) ' 


UK,  Jo) 


k  sin# 


(B5) 


(f  /)"  1  •'  -  — J„UK,n)J.,(k a  sintf ), 

(P . /a)ksu\0 


a>p(e,/p,)''2 


/?„  sin/10 


^(A„.,X/„(kasin«)[l  +  R"(0)]. 


2r  ( Pn,„/o)k  sin# 

Next,  we  shall  evaluate  the  E„  component  which  is  derivable  from  Eqs.  (6)  and  (7).  Thus 

E„  =  ~  ( 1  -  R  " M(# ) )  —  f  f  p  dp' d</>‘ [K,,  cos(0  -  0 ') 

4tt(o£  r  J, i  J ii 

4- A',t  sin(0  -  0 ')]  exp[  —  ikp'  sin#  cos(0  —  0 ')]. 

Letting  the  integral  part  of  the  above  expression  be  I2,  and  substituting  with  (11a)  and  (lib),  we  have 

/,=  f  f  p  dp' d<y  i(e  l//i,)‘''£(Pf;f^-p,)cos/i<^  *  cos(0  -  0 ')  -  ~ — sin/70  ’  sin(0  -0') 
J»  J«  \  a  J  P„„,p  \  a  / 

a  exp|  —  ikp  sin# cos(0  -  0')|. 

It  can  be  shown  that 


r 
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cos/70  '  cos(0  —  0 ')  exp[  -  ikp'  sin#  cos(0  —  0 ')  |  ~  —2 n  (  --/)"' 1  cos/70./,',  (kp'  sin# ), 
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kp  sintV 


Therefore 


■  cos/707,,  (kp  sin# ). 


n'J„(kp'  sin# )J„  [( P„m/a)p' ] 


It  can  be  shown  that 

,,.(P . \,.l,  •  ^  n2Jn[(pn„,/a)p']J„(kp’sin(i) 

p  j p  y -n(kP  sm«)+  - 


(B7) 
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(BIO) 
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As  such 
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The  Elliptical  Microstrip  Antenna  with 
Circular  Polarization 

L.  C.  Shen 

Department  of  Electrical  Engineering 
University  of  Houston 
Houston,  Texas  77004 

Abst  ract 

Theoretical  analysis  of  an  elliptical  microstrip  antenna  is  presented. 
It  is  shown  that  radiation  from  this  antenna  is  left-hand  or  right-hand 
circularly  polarized  in  a  narrow  frequency  band  when  the  eccentricity  of 
the  ellipse  is  small.  The  theoretical  result  agrees  with  the  experimental 
data. 


D-l 


1.  Introduction 


Microstrip  antenna.-;  of  circular  or  rectangular  shapes  usually  radiate 
linearly  polarized  waves.  In  many  practical  applications  circular  polari¬ 
zation  is  required,  circular  polarization  may  he  obtained  by  usinq  multiple 
foods  [11  or  by  altorinq  the  shape  of  a  rectangular  microstrip  antenna  [2). 

It  is  also  known  that  a  simple  microstrip  antenna,  namely,  an  elliptical  patch 
on  a  microwave  printed-circuit  board,  can  be  made  to  radiate  circularly  polar¬ 
ized  waves  (31.  Such  an  antenna  requires  only  one  feed  and  its  geometrical 
shape  is  simple  enough  to  permit  rigorous  theoretical  analysis  to  be  carried 
out  in  a  standard  coordinate  system.  In  the  present  study,  the  radiated  field 
of  an  elliptical  microstrip  antenna  is  obtained  in  terms  of  tabulated  functions. 
The  effects  of  the  fringe  field  at  the  edge  of  the  elliptical  patch  and  those  of 
the  dielectric  substrate  are  taken  into  account  in  the  calculation. 

2.  The  Elliptic  Coordinates 

An  elliptical  microstrip  antenna  is  shown  in  Fig.  1  in  which  the  parameters 
a,  b,  c  and  d  are  also  defined.  The  elliptical  conducting  patch  may  be  obtained 
by  applying  photolithographic  techniques  to  a  microwave  printed-circuit  board. 
The  antenna  may  be  driven  on  the  edge  at  <f>0  degree  from  the  semi-major  axis  by 
a  coaxial  line  through  the  dielectric  substrate  (as  shown  in  Fig.  1)  or  by  a 
microstrip  line. 

Elliptical  coordinates  of  u  and  v  are  used  in  substitution  of  x-y  coordi¬ 
nates  while  the  z  axis  remains  unaffected.  For  an  ellipse  of  semi-major  axis 


a  and  semi-minor  axis  b,  the  foci  are  at  x  =  ±c  where 


2 


c 


(,2-bV/2 


(i) 


The  x-y  coordinates  are  related  to  the  u-v  coordinates  by  the  following  for¬ 
mulas  : 


X  =  C  cosll  11  cosv 

(2a) 

y  =  c  sinh  u  sinv  . 

(2b) 

The  ellipse  is  defined  by 

„  ,a+b, 

u  =  uQ  =  in  (— — )  , 

(3) 

and  the  driving-point  is  at  the  elliptical  coordinates  (u  =  u^,  v  =  v^) 
where 


-1  tan  ^0 
v  =  tan  x(  ■  -  V) 

0  'tanh  u  ' 


(A) 


The  eccentricity  of  the  ellipse  is  defined  as 


(5) 


Assuming  that  the  thickness  of  the  dielectric  substrate  is  much  smaller  than 
a  wavelength  in  the  dielectric  so  that  fields  in  the  dielectric  do  not  vary 
with  respect  to  z,  the  Maxwell's  equations  reduce  to  the  following  set  of 
equations: 


2  2 

9  E  9  E 

- z  +  - z 

2  2 
9  u  9  v 


2  2  2  2 
+  k,c  (cosh  u-cos  v)E 
1  z 


0 


(6) 


H  = 
u 


H  = 
v 


2  2 

coyc/cosh  u-cos  v 


r>E 
_ v 

3v 


,  DE 

ZJ _  _* 

Du 


wpc/cosh^u-cos^v 


H  =  0,  K  =  0,  and  K  -  0, 
?.  u  v 


and 


(7) 

(8) 


ki  = 


(9) 


Throughout  the  analysis  the  time  factor  e^Ut  is  used;  the  magnetic  permea¬ 
bility  p  of  the  dielectric  is  equal  to  that  of  the  air,  that  is,  p  = 
and  the  permittivity  of  the  dielectric  and  that  of  air  are  denoted  and 
C^,  respectively. 

3.  Mathieu  Functions 


Solutions  to  the  partial  differential  Equation  (6)  may  be  expressed  in 


terms  of  Mathieu  functions  S  (C.cosv)  and  S  (C,cosv)  and  related  radial 

om  em 

functions  J  (C,coshu)  and  J  (C,coshu)  [4], 
om  em 


00 

E  =  y  [A  S  (C,cosv)J  (C.coshu) 
z  ‘  n  m  om  om 

m=u 


+  B  S  (C,cosv)J  (C,coshu)l  (10) 

m  em  em 


where  C  =  k^c  . 


In  practical  applications,  the  excitation  of  the  m  =  1  mode  is  favored  since 


4 


it  has  the  lowest  resonant  frequency  and  its  radiation  is  mainly  in  the  z 
direction.  Therefore  and  from  hereafter,  only  two  terms  in  (10)  correspond¬ 
ing  to  the  m  =  1  mode  will  be  retained. 

Since  the  electric  field  has  a  z  component  only,  the  boundary  condition  on 
the  upper  conducting  patch  and  on  the  lower  conducting  ground  plate  is  auto¬ 
matically  satisfied.  The  only  boundary  condition  to  be  imposed  is  on  the 
elliptical  edge  defined  by  u  =  u^: 

-»v  =  Kz+Vs  (at  u  =  Uq>  *  (ll) 


In  the  above  equation,  H  can  be  calculated  from  (8)  and  (10);  K  is  given 

v  7. 

by  (10);  is  the  current  maintained  by  the  transmission  line  at  the  driv¬ 
ing  point  (Ug,vQ) : 


K 

z 


2  2 
C/  cosh  Uq-cos v 


<$(v-vQ) 

0 


(12) 


where  1  is  the  total  current  at  the  driving  point,  and  y  is  the  surface 
u  s 

admittance  as  seen  by  the  elliptical  cavity  bounded  by  u<u^  and  -d<z<0. 
The  surface  admittance  takes  into  account  the  effect  of  radiation  and  the 
fringe  field.  Its  value  will  be  given  later. 

Multiplying  S0^  on  both  sides  of  (11)  and  integrating  over  v  from  0 
to  2tt  yields  the  following: 


A1  =  9J 


-j<4i  I0Sol(C,cos  vQ)/^ 


ol 

du 


2  2 

u  +jygwyC/cosh  u0-cos  vq  Jol(C,cosh  uQ) 


(13) 


may  be  obtained  by  a  similar  procedure: 
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B,  "  ;ij 


-jojM  i0Si,((;,«;msv0)/M| 


(14) 


---I  +  jy  wpC/cosli^u  -cos^v.  J  ,(C,coshu.) 

-)u  'uQ  JJs  0  0  el  0 


where  = 


2tt  2  f2v  „ 

S  .  dv  and  Mf  =  S  dv  . 

o  o1  1  Jo  el 


4.  Approximate  Formulas 


Experimental  data  [3]  showed  that  the  elliptical  antenna  would  produce 
circular  polarization  only  if  the  eccentricity  was  small.  Difference  between 
a  and  b  was  of  the  order  of  a  lew  percent.  That  is,  the  ellipse  should  al¬ 
most  be  a  circle.  For  small  eccentricity,  C  =  k^c^l,  approximate  formulas 
for  the  Mathieu  functions  are  available  [4]: 


Sol(C,cosv)  =  (l+3C2/32)sinv-(C2/32)sin3v 

J  . (C,coshu)  =  ^ - 4 -  [ (l+3C2/32) J  (C  sinhu) 

1  1  (l+CT/8)  1 

-  (C2/32)J3(C  sinhu)] 

S  ,(C,cosv)  =  (1+C2/32)cosv-(C2/32)cos3v 
el 

Jel(C,coshu)  =  /|  ((1+C2/32)J1(C  coshu) 

+  (C^/32)J3(C  coshu)] 

M*  =  7T  ( 1  +  C2/16) 

M°  «  tt(1  +  3C2/16) 


(15) 

(16) 

(17) 

(18) 

(19) 

(20) 


In  (16)  and  (18),  and  are  Bessel  functions  of  order  1  and  3,  respec¬ 
tively.  Again,  since  the  ellipse  is  almost  a  circle  with  0.97  <  b/a  <  1 


the  admittance  surface  y  of  a  circle 

s 


be  used  | 5]: 


ys  =  8s  +  j  bS  * 


gs  =  2a rQ  [(k0a)  X1  +  V 


1.  =  [J.(k.a  sinO)]  sinO  dO  , 

1  '0  1  u 

TT  2„  , 

COS  [I. (k  a  sinO) \  dO, 


X2  " 


sinO  1  r~0 


-  Jj  (w)/  |C  j.J  |  (w)  ]  , 


w  =  1.841/(]+A)l/2  , 

a  E1 

A  =  IMf-r)  +  1.41  +  1.77 

TTC^cl  zd 


+■  -  (0.268  —  +  1.65) ] 
a  f-0 


^0 


1/2 


h  =  W 


1/2 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 
(29) 


The  surface  conductance  gg  is  determined  by  the  power  radiated  through  the 
opening  of  the  cavity,  and  the  surface  susceptance  bg  is  determined  by  the 
fringe  field  which  increases  the  capacitance  by  a  factor  of  (1+A).  The 
formula  for  A  stated  here  is  due  to  Chew  and  Kong  [6]  who  obtained  this 
formula  which  is  more  accurate  than  the  one  used  in  [5], 

5.  Radiated  Field 

The  internal  fields  of  the  elliptical  microstrip  antenna  are  completely 


determined  in  terms  of  the  approximate  formulas  presented  in  the  preceding 
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section.  The  radiated  electric  field  Er  can  then  be  calculated  from  E 

2 

which  sets  up  an  equivalent  magnetic  current  at  u  =  u^,  0  <  v  <  2n: 

Er  =  Vx?  CIO) 

f2v 

i"  -  -  ti  dv' 

^0  v 

2  2 

where  h  =  C/cosh  u  -cos  and  G  is  the  f ree-space  Green's  function.  By 
v  0  0 

using  the  f ree-space  Green's  function  the  effect  on  the  radiation  field  of 
the  dielectric  layer  in  the  external  region  (u  >  u^)  is  neglected.  However, 
according  to  an  analysis  by  Kong  [7],  the  effect  can  be  taken  into  account 
by  introducing  a  factor  to  the  final  result.  This  will  be  done  later  in 
the  present  analysis. 

Substituting  (10)  in  (31),  noting  that 

A  ^  A 

h  v'  =  -a  sinv  x  +  b  cosv  y  ,  (32) 

v 


and  after  rather  tedious  manipulations  of  the  integral,  the  following  re¬ 
sults  are  derived: 

"jkor  k2 

Eo  =  -  A  -  --  cos9[l-R™(0)  ]}(I  sin<t>-I  cos$) 

0  47tr  k2-k2sin20  X  y  (33) 

“jkor 

E.  =  -  [l+RTE(e)](I  cos0+I  sin<J>)  (34) 

$  4rrr  x  y 

where  k^  and  k^  are  wavenumbers  In  free  space  and  in  the  dielectric,  re¬ 
spectively.  Spherical  coordinates  r,  0,  <j>  are  used.  Furthermore, 


!-K™(0) 


-,2klxd 

)-o _ ,,_R™ 

■  TM  -  |2k  i  xd  U  or  ’ 
1+RI)1C 


i+rte(0) 


->2klzd 

l-e  _  n+„TE) 

01  ' 

1  K0lt- 


k.  —  lc 

TM  0 z  lz 

01  '  kOz+kU  ’ 


TE  _  k0.z_'lC-'<iJ— ~- 

*01  k  k’f+k?k 

()z  1  0  lz 


k0z  ”  k0COS°  * 


k,2  -  (k:’-k„  Sln20)*/2 


I  =  --nafA1  [JQ-J2c°s2a  +  ( 3J0-4J2cos2a+Jj!|CosAa)  ] 


p  ^ 

-  B'[J2sin2a  -  JAsin4a]} 


I  =  TTb{A'  [-J?sin2a  +  ^  (_2J2sin2a+J^sin4a) ] 


C  " 

+  B’[J0+J2cos2a  +  32  (JQ-J4cos4a) ] } 


A*  =  Ax  Jq1(c»  cosh  uq)  • 


B’  *  Bl  Jel(c»  cosh  u0>  * 


-1  ,a  cos4>s 

°  *  tan  (rsiS)’ 


and  J  ,  J  ,  and  J.  are  Bessel  functions  of  under  0,  2,  and  4,  respectively 
0  2  4 


with  argument  <)  where 


22  22  1/2 
Q  =  k.QSinO(a  cos  $  +  b  sin  (j>) 


(46) 


In  (33-34)  factors  have  been  introduced  to  take  into  account  the  presence 
of  the  dielectric  layer. 


6.  Results  and  Discussion 

The  radiated  electric  fields  are  qiven  in  (33)  and  (3d)  are  in  terms  of 
simple  trigonometric  functions  and  Bessel  functions  for  which  efficient  com¬ 
puter  programs  are  readily  available.  In  Fig.  2  the  amplitude  and  phase  of 
E^  and  on  the  z-axis (which  is  the  direction  of  maximum  radiation)  are 
plotted  as  functions  of  frequency,  with  the  position  of  the  driving  point 
as  a  variable  parameter.  It  is  seen  that  optimum  circular  polarization 
can  be  obtained  when  the  driving  point  is  placed  at  43°  to  the  major  axis 
with  frequency  equal  to  2.79  GHz.  The  effect  of  changing  the  position  of 
the  driving  point  is  to  change  the  relative  amplitudes  of  the  x  directed 
and  the  y-direction  electric  fields.  Moving  the  driving  point  toward  the 
major  axis  enhances  the  excitation  of  current  flowing  parallel  to  the  major 
axis  and  thus  results  in  greater  amplitude  for  E^.  It  is  also  seen  that 
the  phase  is  not  a  sensitive  function  of  4>q.  The  elliptical  microstrip 
antenna  for  which  Fig.  2  is  produced  has  b/a  =  0.98.  In  Fig.  3,  the  fields 
of  a  microstrip  antenna  of  b/a  =  0.96  are  calculated.  It  is  seen  that 

although  E  is  equal  to  E  in  amplitude  at  f  =  2.82  GHz,  the  phase  differ- 
x  y 

t 

ence  is  greater  than  90°  so  that  the  quality  of  circular  polarization  is 
not  as  good  as  that  shown  in  Fig.  2.  Note  that  in  Fig.  2  as  the  frequency 


is  chanqed  from  the  optimum  frequency,  not  only  do  the  relative  ampliudes 
of  Ex  and  Ey  chanqe  but  also  the  relative  phase  angles  deviate  from  the 
ideal  90°  value.  Polarization  ellipses  on  the  z  axis  are  shown  in  Fiq.  4 
as  functions  of  the  operating  frequency.  Arrows  indicate  direction  of 
rotation  of  !•:  field  at  a  fixed  point  in  t:h<-  z  axis.  Thus  in  Fie.  4  all 
fields  are  left-hand  el  1  i  pt  i<~,i  1  ly  polar  i  zed . 

Radiation  patterns  oi  an  elliptical  antenna  with  b/a  =  0.9H  operated 
at  the  optimum  frequency  for  circular  polarization  are  shown  in  Figs.  5a 
and  5b.  It  is  seen  that  quality  of  circular  polarized  is  quite  good  for  0 
less  than  45°.  The  polarization  ellipses  are  shown  in  Fig.  6  for  various 
6  angles. 

Based  on  extensive  experimental  data,  Yu  [3]  discovered  the  following 
characteristics  of  the  elliptical  microstrip  antenna. 

(a)  "Generally,  the  radiation  associated  with  the  elliptical  antenna 
element  is  elliptically  polarized,  but  is  circularly  polarized 
when  the  antenna  element  is  coupled  through  a  feed  point  on  a 
radial  line  of  the  elliptical  lamina  which  is  or<ented  at  a  45° 
azimuthal  angle  relative  to  a  semimajor  axis  of  the  ellipse.  If 
the  radial  line  on  which  the  feed  point  is  located  is  rotated  45° 
counter  clockwise  relative  to  a  semimajor  axis  of  the  ellipse, 
the  polarization  of  the  radiation  communicable  by  the  antenna 
element  will  be  lef t-circularly  polarized.  If  the  radial  line 
on  which  the  feed  point  is  located  is  oriented  at  an  angle  of  45° 
clockwise  relative  to  a  sem.imajor  axis  of  the  ellipse,  the  polar¬ 
ization  of  the  radiation  in  question  will  be  r ight-circularly 
polarized".  These  statements  are  in  agreement  with  the  present 
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theory,  as  shown  in  Figs.  4  and  6.  1 1  is  also  noted  that  if  0^ 

is  changed  to  it-0^,  and  0  to  ti-0  in  all  preceding  formulas,  then 

E  will  be  unchanged  while  E,  is  changed  to  -F. ,  ,  resulting  in 
u  99 

right-hand  clliptically  polarized  radiation. 


(b)  To  achieve  an  operating  frequency  f,  the  semimajor  axis  a  of  the 
ellipse  is  given  by  the  equation  a  -  p/f/pt;^,  where  p  is  an 
empirical  constant  ranging  from  0.27  to  0.29.  Accordi  lg  to  the 
present  theory,  using  parameters  in  Figs.  2,  5,  and  6,  p  is  found 
to  be  equal  to  0.275  which  agrees  very  well  with  the  empirical 
value. 


(c)  "It  has  been  found  that  the  desired  circularly  polarized  radiation 
communicability  may  best  be  achieved  by  limiting  the  eccentricity 
of  the  ellipse  to  a  range  of  10  to  20%. "  For  b/a  =  1.84/1.88, 
the  eccentricity  is  equal  to  20.5%.  Present  theoretical  results 
show  that  good  circular  polarization  is  obtained  [Fig.  4(a)].  For 
b/a  =  1.80/1.88,  the  eccentricity  is  equal  to  28.9%.  Present 
theory  shows  that  quality  of  circular  polarization  is  not  very 
good,  as  seen  from  Fig.  4(b). 

7.  Conclusion 


Currents,  internal  fields  and  radiated  fields  of  an  elliptical  micro¬ 
strip  antenna  have  been  expressed  in  terms  of  Mathieu  functions  which  may 
be  expressed  in  terms  of  Bessel  functions  when  the  eccentricity  of  the 
ellipse  is  small.  Left-hand  or  right-hand  circular  polarization  can  be 
obtained  in  a  narrow  frequency  band  for  an  elliptical  microstrip  antenna 


of  small  eccentricity. Only  one  feed  is  needed  to  excite  the  circularly  polar 
ized  antenna.  All  formulas  are  in  terms  of  functions  that  can  easily  be 
programmed  on  small  computers  to  facilitate  the  design  of  such  an  antenna 
for  practical  applications.  Theoretical  results  are  in  good  agreement  with 
available  experimental  data. 
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Figure  Legends 


Fig.  1  An  elliptical  mlcrn.st  rip  antenna  and  the  elliptical  coordinates. 

Fig.  2  Amplitudes  and  phases  of  F.  and  E  on  the  z  axis  of  an  elliptical 

x  y 

microstrip  antenna  with  a  =  1.88  cm,  b  =  1.84  cm,  d  =  0.16  cm, 

^  =  2.47  Cq,  driven  at  with  1^  =lamp. 

Fig.  3  Amplitudes  and  phases  of  E  and  E  on  the  z  axis  of  an  elliptical 

x  y 

microstrip  antenna  with  a  =  1.88  cm,  b  =  1.80  cm,  d  =  0.16  cm, 
tj  =  2.47  Lq,  driven  at  4>q  with  1^=1  amp. 

Fig.  4  Polarization  ellipses  on  the  z  axis  as  functions  of  frequency 

(a)  for  elliptical  antenna  with  b/a  =  1.84/1.88  (b)  for  elliptical 

antenna  with  b/a  =  1.80/1.88.  Other  parameters  are  same  as  those 
given  in  the  legend  of  Fig.  2.  Arrows  indicate  direction  of  rota¬ 
tion  of  the  electric  field.  Antennas  are  driven  at  ■  45°, 
resulting  in  left-hand  circular  polarization. 

Fig.  5  Radiation  patterns  of  an  elliptical  microstrip  antenna  with  para¬ 
meters  given  in  the  legend  of  Fig.  2.  Driving-point  is  at  $q  “  45° 

(a)  0  =  0  and  it  plane  (b)  0  =  tt/2  and  3tt/2  plane  f  =  2.79  GHz. 

Polarization  ellipses  as  functions  of  9  produced  by  an  elliptical 
microstrip  antenna  with  parameters  given  in  the  legend  of  Fig.  2, 
driven  at  <p  “  45°,  with  f  =  2.79  GHz  (a)  on  the  4>  *  0  plane 

(b)  on  the  <f>  =  90°  plane. 
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Abstract 

Elliptical  shaped  printed-circuit  antennas  were  fabricated  with 
varying  eccentricities  and  their  impedance  and  radiation  patterns  were 
measured.  Special  attention  was  devoted  to  investigating  the  design 
criteria  which  produced  the  best  circular  polarization. 
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I.  Introduction 


Over  the  past  several  years  the  printed-circuit  antenna  has  proven  to 
he  an  efficient  and  effective  radiator.  In  its  simplest  configuration, 
whether  rectangular  or  circular,  it  can  he  designed  to  produce  a  broad 
beamwidth,  linearly  polarized  radiation  pattern  with  its  maximum  in  the 
direction  normal  to  the  plane  of  the  antenna.  The  most  direct  approach 
to  provide  circular  polarization  from  such  an  antenna  is  to  use  two  feeds 
located  geometrically  90  degrees  apart  and  with  a  relative  phase  shift  of 
90°.  This  configuration  excites  two  orthogonal  modes  each  of  which  pro¬ 
vides  a  linearly  polarized  wave  at  right  angles  to  each  other  and  shifted 
in  phase  by  90°. 

Several  methods  have  been  proposed  to  provide  the  desired  circular 
polarization  without  the  additional  complexities  inherent  in  the  dual  feed 
and  phase-shifter  radiator.  The  pentagon  [1],  corner-red  rectangle  [2], 
and  the  45°  slot  [JJ  variations  have  all  been  designed  for  this  purpose. 
Recently  Shen  [4]  has  shown  that  circular  polarization  could  also  be 
expected  from  a  slightly  elliptical  radiator  fed  along  a  line  45°  from  its 
major  axis.  Similar  behavior  had  been  seen  experimentally  by  Yu  [5]  earlier 
but  never  developed. 

In  this  paper,  a  careful  experimental  study  is  carried  out  to  investi¬ 
gate  the  behavior  of  the  radiated  fields  and  the  impedance  of  the  ellipti¬ 
cal  radiator  as  a  function  of  frequency,  eccentricity,  and  dielectric  thick¬ 
ness.  The  results  are  then  compared  with  the  theoretical  predictions  of 
Shen  [4J. 
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II.  Kxper imcntal  Impedance  Measurements 

To  Investigate  the  properties  of  the  elliptical  printed-circuit  antenna 
a  set  of  ellipses  with  varying  eccentricities  were  etched  on  two  different 
thicknesses  of  printed-circuit  board.  The  slightly  elliptical  discs  shown 
in  Figure  1  varied  from  a  b/a  ratio  of  1.00  (that  of  a  circle)  to  b/a  *=  0.96. 

The  major  radius  remained  constant  at  4  cm  while  the  minor  radius  was  pro¬ 

gressively  reduced  to  about  3.84  cm.  Five  different  ellipses  were  etched  on 
both  0.3175  cm  (~l/8")  and  0.1575  cm  (~1/16")  thick  teflon-fiberglass  micro- 
wave  printed-circuit  board. 

The  impedance  of  each  antenna  was  measured  as  a  function  of  frequency 
from  1  to  2  GHz  using  an  automatic  network  analyzer.  An  illustrative  example 

of  one  of  these  swept  frequency  measurements  is  shown  in  Figure  2  in  Smith 

chart  format  for  the  case  of  d  =  .3175  cm  andb/a=.976.  While  no  attempt 
was  made  to  match  the  antennas  to  the  characteristic  impedance  of  the  50  ohm 
line,  this  could  be  accomplished  by  feeding  the  disc  away  from  the  edge 
along  the  same  45°  radial  line.  The  experimental  results  are  also  shown  in 
a  display  of  the  real  and  imaginary  parts  of  the  input  impedance  as  a 
function  of  frequency.  The  resistance  and  reactance  are  shown  for  the  case 
of  d  =  .3175  cm  for  four  different  values  of  b/a  in  Figurts  3-b.  Similar 
graphs  for  the  case  of  d  =  .1575  cm  are  shown  in  Figures  7-10. 

The  deviation  from  a  circle  seen  in  Figure  2  is  manifest  in  the  double 
peaking  behavior  seen  for  the  same  antenna  in  Figure  5.  The  impedance  of 
the  radiators  that  are  most  nearly  circular  (Figures  3  and  7)  are  seen  to 
resemble  the  impedance  of  a  circular  disc  as  one  would  expect.  As  the  disc 
becomes  more  elliptical  (Figures  4  and  8)  two  distinct  resonances  begin  to 
appear.  This  feature  is  more  apparent  in  the  0.1575  cm  thickness  case 


since  tiie  resonance  curve  for  the  circular  disc  of  this  thickness  is  approxi¬ 
mately  half  as  wide  as  for  the  thicker  case.  Even  though  the  resonance 
curves  in  Figures  3-10  are  seen  to  be  broader  in  the  double  resonance  cases 
the  polarization  does  not  remain  circular  over  this  larger  impedance  band¬ 
width.  At  the  first  of  the  double  peaks  the  radiation  is  linearly  polarized 
along  the  major  axis  of  the  ellipse.  As  the  frequency  is  increased  the 
radiation  becomes  more  nearly  circularly  polarized.  As  the  second  peak  is 
approached  the  radiation  becomes  linearly  polarized  again,  but  this  time 
oriented  along  the  minor  axis  of  the  ellipse. 

III.  Radiation  Patterns 

To  investigate  the  polarization  of  the  radiated  fields  from  these  ellip¬ 
tical  printed-circuit  antennas,  far  field  radiation  patterns  were  measured. 

Since  the  possibility  of  circular  polarization  was  of  concern,  all  patterns 
were  taken  with  a  spinning,  linear-polarized  horn.  The  printed-circuit  an¬ 
tenna  itself  was  used  as  the  source  by  attaching  a  small,  lightweight, 
battery-powered  transmitter  directly  to  the  back  of  the  ground  plane  of  the 
antenna.  This  transmitter  assembly  was  placed  on  the  rotating  table  inside 
an  anechoic  chamber  and  the  spinning  horn  then  used  as  the  receiver. 

Radiation  patterns  are  shown  in  Figures  11-13  for  three  different  ratios 
of  b/a  with  the  same  thickness  d  =  .3175  cm.  The  4>  =  90°  plane  represents 
a  cut  directly  over  the  minor  axis  of  the  elliyce.  Similar  patterns  were 
also  recorded  for  cuts  made  over  the  major  axis  and  over  the  feed  point  of 
the  disc.  The  frequencies  were  chosen  by  experimentally  finding  the  fre¬ 
quency  for  which  the  axial  ratio  was  a  minimum  at  0  =  0°  (normal  to  the 
plane  of  the  disc). 
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For  each  case  the  spin-linear  measured  patterns  are  shown  by  the  solid 

line.  In  addition  the  theoretical  maximum  and  minimum  fields  for  all  pos¬ 

sible  orientations  of  the  linear  polarized  receiving  antenna  are  shown  for  this  fre¬ 
quency  and  geometry.  The  curve  is  normalized  only  at  the  point  of  the 
maximum  field  at  0  =  0.  It  is  seen  that  both  the  theoretical  curves  and 
the  experimental  data  confirm  that  for  a  properly  chosen  eccentricity  and 
frequency  the  resulting  radiation  is  quite  nearly  circularly  polarized  (about 
1  dB  for  b/a  =  .976  in  Figure  12).  This  circularity  is  seen  to  deteriorate 

both  for  the  case  when  the  disc  is  made  more  circular  (Figure  11)  and  when 

it  is  made  more  elliptical  (Figure  13). 

To  investigate  its  dependence  on  frequency  the  axial  ratio  at  0  =  0° 
was  measured  for  each  antenna  as  the  frequency  was  varied.  Using  the  exact 
dimensions  of  the  experimental  antennas  and  the  manufacturer's  stated  value 
of  dielectric  constant  (Cf  =  2.48)  the  theoretical  curves  accurately  indicate 
the  lowest  value  of  axial  ratio  and  the  general  shape  of  the  curves  for  each 
eccentricity.  However,  the  frequency  for  which  the  minimum  axial  ratio 
occurs  is  not  well  predicted.  For  this  reason  a  search  for  the  indicated 
value  of  dielectric  constant  was  made  using  the  theoretical  formulas  and 
the  experimentally  found  frequency  of  minimum  axial  ratio.  For  each  differ¬ 
ent  antenna,  the  best  value  was  found  to  be  e  =  2.41.  Thus  in  Figure  14, 
the  experimental  data  and  the  theoretical  curves  for  this  value  of  dielectric 
constant  are  shown.  Reasonable  correlation  is  seen  for  all  aspects  of  the 
theory  and  experiment.  (Use  of  =  2.48  would  simply  shift  the  theoretical 
curves  to  the  right  along  the  frequency  axis.) 


These  curves  indeed  show  that  there  exists  an  optimum  eccentricity  and 
also  give  an  indication  of  possible  frequency  bandwidth.  As  noted  earlier. 
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this  bandwidth  due  to  polarization  is  even  narrower  than  the  impedance  band¬ 
width  measured  in  the  previous  section.  Kven  for  the  best  antenna 
(b/a  =  .976),  the  axial  ratio  is  less  Lhnn  6  dB  only  over  the  range  from 
L.33  to  1.36  GHz  or  a  bandwidth  of  1.5%. 

IV.  Conclusions 

It  has  been  shown  that  with  the  proper  choice  of  eccentricity  and  fre¬ 
quency  the  elliptical  printed-circuit  antenna  can  provide  a  circularly 
polarized  radiation  field  using  only  a  simple  single  feed  network.  The 
behavior  of  its  impedance,  pattern  and  axial  ratio  as  a  function  of  eccen¬ 
tricity,  frequency,  and  dielectric  thickness  has  also  been  investigated. 

It  is  seemingly  an  excellent  design  choice  when  both  the  usual  character¬ 
istics  of  printed-circuit  antennas  and  circular  polarization  are  required. 
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Figure  1.  Geometry  of  elliptical  printed-circuit  antenna 

Figure  2.  Smith  chart  display  of  impedance  (d  =  .3175  cm;  b/a  *  . 

Figure  3.  Impedance  versus  frequency  (d  =  .3175  cm;  b/a  =  .993) 

Figure  4.  Impedance  versus  frequency  (d  =  .3175  cm;  b/a  *  .985) 

Figure  5.  Impedance  versus  frequency  (d  =  .3175  cm;  b/a  =  .976) 

Figure  6.  Impedance  versus  frequency  (d  =  .3175  cm;  b/a  =  .960) 

Figure  7.  Impedance  versus  frequency  (d  =  .1575  cm;  b/a  =  .996) 

Figure  8.  Impedance  versus  frequency  (d  =  .1575  cm;  b/a  -  .983) 

Figure  9.  Impedance  versus  frequency  (d  =  .1575  cm;  b/a  =  .976) 

Figure  10.  Impedance  versus  frequency  (d  *  .1575  cm;  b/a  ■  .962) 

Figure  11.  Radiation  pattern  (d  =  .3175  cm;  b/a  *  .985) 

Figure  12.  Radiation  pattern  (d  *  .3175  cm;  b/a  *  .976) 

Figure  13.  Radiation  pattern  (d  *  .3175  cm;  b/a  =  .960) 

Figure  14.  Axial  ratio  versus  frequency  (d  =  .3175  cm) 
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